The mass ejection from the merger of binary neutron stars 
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Numerical-relativity simulations for the merger of binary neutron stars are performed for a variety 
of equations of state (EOSs) and for a plausible range of the neutron-star mass, focusing primarily 
on the properties of the material ejected from the system. We find that a fraction of the material 
is ejected as a mildly relativist ic and mildly anisotropic outflow with the typical and maximum 
velocities ^ 0.15 - 0.25c and ~ 0.5 - 0.8c (where c is the speed of light), respectively, and that the 
total ejected rest mass is in a wide range 10~^ - lO^^M©, which depends strongly on the EOS, the 
total mass, and the mass ratio. The total kinetic energy ejected is also in a wide range between 
10^^ and 10^^ ergs. The numerical results suggest that for a binary of canonical total mass 2.7 Mq^ 
the outflow could generate an electromagnetic signal observable by the planned telescopes through 
the production of heavy-element unstable nuclei via the r-process pil3j or through the formation 
of blast waves during the interaction with the interstellar matter U, if the EOS and mass of the 
binary are favorable ones. 

PACS numbers: 04. 25. Dm, 04.30.-w, 04.40.Dg 



I. INTRODUCTION 

Coalescence of binary neutron stars is one of the 
most promising sources for next-generation kilo-meter- 
size gravitational- wave detectors such as advanced LIGO, 
advanced VIRGO, and KAGRA (LCGT) [5 . These de- 
tectors will detect gravitational waves in the next 5 - 
10 yrs. Statistical studies have predicted that the de- 
tection rate of gravitational waves emitted by binary 
neutron stars for these detectors will be ^ 1 - 100 per 
year [6l [7]. The typical signal-to- noise ratio for most 
of these events will be 10 or less. Thus, it will be 
quite helpful if electromagnetic or other signals observ- 
able are associated with the gravitational-wave bursts 
and the gravitational-wave detection is accompanied by 
the detections of other signals. Short-hard gamma-ray 
bursts (SGRB) have been inferred to be associated with 
the binary neutron star merger IQ. However, the jet 
of SGRB would be highly collimated [9|, and hence, it 
will not be always possible to detect SGRB associated 
with the binary neutron star mergers. Moreover, it is 
not guaranteed that the telescopes for the observation of 
SGRB will be in operation with the gravitational- wave 
detectors. Exploring other possible signals that could 
be detected is a very important subject in the fields of 
gravitational- wave physics/astronomy [3l |4l IToUTSl 118] . 

This paper presents our latest results of numerical sim- 
ulations performed in the framework of numerical rela- 
tivity, focusing in particular on the exploration of the 
material ejected from binary neutron star mergers. In 
the past decade, numerical simulations for the merger of 
binary neutron stars in full general relativity, which is 
probably the unique approach of the rigorous theoretical 



study for this subject, have been extensively performed 
since the first success in 2000 [21 (see, e.g., [22j [23] for 
a review of this field). However, most of the simulations 
have focused on the studies of gravitational waveforms 
and the resulting product formed in the central region. 
Few attention has been paid to the study for the mate- 
rial ejected (but see ^3 for a study in an approximate 
framework of general relativity, and see p!4UT7l for an 
early effort in the context of Newtonian gravity) . 

The material ejected from binary neutron star merg- 
ers may generate electromagnetic signals observable in 
the current and future-planned telescopes. One possi- 
ble signal could be generated by the radioactive decay of 
unstable r-process nuclei, which are produced from the 
neutron-rich material in the ejecta [JHS [EJ |13l HT] [T8] . 
A fraction of the unstable nuclei produced subsequently 
decay in a short timescale and could heat up the ejecta, 
which emits a UV and visible light that may be observ- 
able by current and future-planned optical telescopes. In 
this case, the typical duration of a peak luminosity is ex- 
pected to be several hours to a day. Another possible 
signal could be generated during the free expansion and 
the subsequent Sedov phase of the ejecta which sweeps 
up the interstellar medium and forms blast waves |4j. In 
this process turning on, the shocked material at the blast 
waves could generate magnetic fields and accelerate par- 
ticles that emit synchrotron radiation in the radio-wave 
band, for a hypothetical amplification of the electromag- 
netic field and a hypothetical electron injection. It is 
also pointed out that the binary neutron star merger 
could drive ultra-relativistic outflows in every direction 
and emit synchrotron radiation in x-ray-to-radio bands 
within a second-to-day timescale [19 . All these stud- 
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ies illustrate that exploring the process of the material 
ejection from binary neutron star mergers in detail is an 
important subject. 

For the detailed numerical study of the ejected mate- 
rial, we have to be careful when following the motion of 
the materials in a low-density outer region. Most of the 
numerical-relativity simulations of binary neutron star 
mergers so far have been performed with a computa- 
tional domain that was not wide enough for this pur- 
pose [22[ . We have to enlarge the computational do- 
main sufficiently widely to conffim that the outflowed 
material is indeed ejected from the system (i.e., we have 
to conffim that the material is indeed unbound by the 
system by following the motion of the ejected material 
for a long time). Another subtle issue in the hydrody- 
namics simulations is that we have to put an artificial 
atmosphere when employing a conservative shock cap- 
turing scheme that is a standard one in this field [24] . In 
our previous simulations [25H27] . we put an atmosphere 
with fairly large density 10^ g/cm^) that did not af- 
fect the motion of neutron stars but did for the motion 
of the ejected material of low density which might es- 
cape to a far region. For the study of the mass ejection, 
we have to reduce the density of the atmosphere as low 
as possible (which should be much lower than the den- 
sity of the ejected material), and in addition, we have to 
carefully conffim that such an artificial atmosphere does 
not affect the properties of the ejected material. In the 
simulation reported in this paper, we succeed in the sim- 
ulation reducing the atmosphere density to a low level 
10^ g/cm^) enough to obtain a scientifically quanti- 
tative result. 

The paper is organized as follows: In Sec. II, we 
summarize the equations of state (EOSs) employed and 
models of binary neutron stars. In Sec. Ill, we briefly 
summarize our formulation and numerics for solving Ein- 
stein's equation and hydrodynamics equations as well as 
the tools for diagnostics. In Sec. IV, numerical results 
are presented, focusing on the properties of the material 
ejected from the system. Section V is devoted to a sum- 
mary and discussion. Throughout this paper, we employ 
the geometrical units c = 1 = G where c and G are the 
speed of light and gravitational constant, respectively, 
although we recover c when we need to clarify the units. 



II. EQUATIONS OF STATE AND CHOSEN 
MODELS 

In this section, we summarize the model EOSs em- 
ployed in this paper, and initial condition of binary neu- 
tron stars chosen for numerical simulations. As shown 
in Sec. \P7\ the properties of the material ejected from 
binary neutron star mergers depend strongly on these 
inputs. 



A. Equations of state 

The exact EOS for the high-density nuclear matter is 
still unknown [28] . This implies that a numerical simula- 
tion employing a single particular EOS, which might not 
be correct, would not yield a scientific result. A study, 
systematically employing a wide possible range of EOSs, 
is required for binary neutron star mergers. Neverthe- 
less, the latest discovery of a high-mass neutron star PSR 
J1614-2230 with mass 1.97 ± O.O4M0 [29 significantly 
constrains the model EOS to be chosen, because it sug- 
gests that the maximum mass for spherical neutron stars 
for a given EOS has to be larger than ~ 2Mq. This indi- 
cates that the EOS should be rather stiff, although there 
are still many candidate EOSs. 

To model a variety of the candidate EOSs, specifically, 
we employ a piecewise poly tropic EOS proposed by Read 
et al. [30 . This EOS is described assuming that neu- 
tron stars are cold (in a zero-temperature state), i.e., the 
rest-mass density, p, determines all other thermodynam- 
ical quantities. To systematically model nuclear-theory- 
based EOSs at high density with a small number of pa- 
rameters, the pressure is written in a parameterized form 
as 

P{p) = nip^' for Pi < p < Pi+i (0<z<n), (1) 

where n is the number of the pieces used to parameterize 
a high-density EOS, pi is the rest-mass density at the 
boundary of two neighboring (i — l)-th and z-th pieces, 
Ki is the polytropic constant for the i-th piece, and 
is the adiabatic index for the i-th piece. Here, po = 0, 
pi denotes a nuclear density ~ 10^^ g/cm^ determined 
below, and Pn+i oo. Other parameters (pi^hCi^Ti) are 
determined by fitting with a nuclear-theory-based EOS. 
Requiring the continuity of the pressure at each p^, 2n 
free parameters, say {f<ii^Ti)^ determine the EOS com- 
pletely. The specific internal energy, e, and hence the 
specific enthalpy, /i, are determined by the first law of 
thermodynamics and the continuity of each variable at 
boundary densities, pi. 

Read et al. [30] showed that a piecewise polytropic EOS 
with three pieces above the nuclear density (i.e., n = 3) 
approximately reproduces most properties of the nuclear- 
theory-based EOS at high density, and they derived the 
ffited parameters for a large number of nuclear-theory- 
based EOSs. In this paper, thus, we employ this piece- 
wise polytropic EOS, determining the free parameters 
basically following j3TH33] (in which a piecewise poly- 
trope with n = 1 was used). First, the EOS below the 
nuclear density pi is fixed by the following parameters 

Fo = 1.35692395, (2) 
K.o/c^ = 3.99873692 x 10"^ (g/cm^)^-^^. (3) 

The EOS for the nuclear matter was determined in [30] 
as follows: p2 was fixed to be p2 = 10^^"^g/cm^, and P2 
at p = p2 was chosen as a free parameter. The reason 
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TABLE I: Parameters and key quantities for four piecewise polytropic EOSs employed in this paper. P2 is shown in units of 
dyn/cm^. Mmax is the maximum mass along the sequences of spherical neutron stars (cf. Fig. [2]). (i?i.35, pi.35) and (i?i.5,pi.5) 
are the circumferential radius and the central density of 1.35M0 and l.SM© neutron stars, respectively. We note that the 
values of the mass, radius, and density listed are slightly different from those obtained in the original tabulated EOSs (see the 
text for the reason). MSI is referred to as this name in |30], but in other references (e.g., [28]), it is referred to as MSO. We 
follow [30] in this paper. 



EOS (log(P2),ri,r2,r3) M^ax(M0) i^i.35(km) pi.35(g/cm^) i^i.5(km) pi.5(g/cm^) 



APR4 (34.269,2.830,3.445,3.348) 


2.20 


11.1 


8.9 X 10'^ 


11.1 


9.6 X 10'^ 


ALF2 (34.616,4.070,2.411,1.890) 


1.99 


12.4 


6.4 X 10^^ 


12.4 


7.2 X 10^^ 


H4 (34.669, 2.909, 2.246, 2.144) 


2.03 


13.6 


5.5 X 10^^ 


13.5 


6.3 X 10^^ 


MSI (34.858,3.224,3.033,1.325) 


2.77 


14.4 


4.2 X 10^^ 


14.5 


4.5 X 10^^ 



is that P2 is closely related to the radius and deforma- 
bility of neutron stars [34 . Namely, P2 primarily deter- 
mines the stiffness of an EOS. Second, pa was fixed to 
be p3 = 10^^-^g/cm^. With these choices, the set of free 
parameters becomes (P2, Fi, Fs). These four param- 
eters are determined by a fitting procedure (see [30j for 
the fitting procedure). 

With the given values of Fi and P2, ni and pi are 
subsequently determined by 



1^1 P2P2 \ 

Pi = (A.oM)^/^^-^°^ 

By the same method, i<i2 and tvs are determined from 



(4) 
(5) 

1 

(6) 



Table [T| lists the EOSs and their parameters which we 
employ in this study. We choose four types of the rep- 
resentative EOSs. APR4 was derived by a variational 
method with modern nuclear potentials [35^ for the hy- 
pothetical components composed of neutrons, protons, 
electrons, and muons; MSI was derived by a mean- field 
theory for the hypothetical components composed of neu- 
trons, protons, electrons, and muons, as well [36]; H4 was 
derived by a relativistic mean-field theory including ef- 
fects of hyperons [37 ; ALF2 is a hybrid EOS which de- 
scribes a nuclear matter for a low density and a quark 
matter for a high density with the transition density is 
3pnuc where pnuc ^ 2.8 x 10^^ g/cm^ [38 . We note that 
the piecewise polytropic EOSs are slightly different from 
the original tabulated ones, because of their simple fitting 
formula. This results in a small error in the mass and ra- 
dius of neutron stars. However, as shown in , the error 
is small (at most several percent), and the semiquantita- 
tive properties of the original EOSs are well captured by 
these simple EOSs. 

Figure [l] plots the pressure as a function of the rest- 
mass density for four EOSs. APR4 has relatively small 
pressure for pi < p ^ ps while it has high pressure for p ^ 
ps. Thus, for p < ps^ which neutron stars of canonical 
mass 1.3 - 1.4M0 have, this EOS is soft, and hence, the 
value of P2 is relatively small. We note that for a small 
value of P2, F2 and/or F3 have to be large (~ 3) because 
the maximum mass of spherical neutron stars, M^ax for 
a given EOS has to be larger than ^ 2Mq. Thus, an 
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FIG. 1: Pressure as a function of the rest-mass density for 
APR4, ALF2, H4, and MSI EOSs (the solid, dashed, dotted, 
and dash-dotted curves, respectively). 



EOS that is soft at p = p2 has to be in general stiff for 
p ^ P3. By contrast, H4 and MSI have pressure higher 
than APR4 for p ^ pa, although the EOSs become softer 
for a high-density region p^ pz- In particular, MSI has 
extremely high pressure (i.e., a higher value of P2) among 
many other EOSs for p ^ pa, and thus, it is the stiffest 
EOS as far as the canonical neutron stars are concerned. 
ALF2 has small pressure for p < p2 as in the case of 
APR4, but for P2 ^ P < Ps, the pressure is higher than 
that for APR4. For p > p2 the pressure of ALF2 is as 
high as that for H4. All the properties mentioned above 
are reflected in the radius, i^i.35, and central density, 
Pi. 35, of spherical neutron stars with the canonical mass 
M = I.35M0 where M is the gravitational (Arnowitt- 
Deser-Misner; ADM) mass of the spherical neutron stars 
in isolation: see Table [l[ The pressure at p = p2 (^2) is 
correlated well with this radius and central density (see 
below). 

Here, a word of caution is necessary for our APR4. 
The pressure in this piecewise polytropic EOS is ex- 
tremely (unphysically) high in the high-density region 
with p ^ 10^^ g/cm^. This results pathologically in the 
situation that the sound velocity exceeds the speed of 
light for the high-density state. In reality, such a high 
density is achieved only in the formation of a black hole 
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(i.e., inside the horizon), and such a pathology may not 
affect the evolution of the system for the outside of the 
horizon. However, this pathology could still break a nu- 
merical simulation after the formation of a black hole. 
To avoid this happens, we artificially set the maximum 
density as 10^^ g/cm^ when employing this EOS. 

Figure |2] plots the gravitational mass as a function of 
the central density and as a function of the circumfer- 
ential radius for spherical neutron stars for four EOSs. 
All the EOSs chosen are stiff enough that the maxi- 
mum mass is larger than 1.97 Mq. Because the pres- 
sure in a density region p ^ 10^^ g/cm^ is relatively 
small (i.e., P2 is small) for APR4 and ALF2, the ra- 
dius for these EOSs is relatively small as ~ 11 km and 
12.5 km, respectively, for the canonical mass of neutron 
stars 1.3 - 1.4M0 [40,. By contrast, for H4 and MSI 
for which P2 is relatively large, the radius becomes a 
relatively large value 13.5 - 14.5 km for the canonical 
mass. The radius has also the correlation with the cen- 
tral density p^. For APR4 and ALF2 with M = 1.35M0, 
Pc ^ 8.9 X 10^"^ g/cm^ and p^ ^ 6.4 x 10^"^ g/cm^. For H4 
and MSI with M = 1.35M0, the central density is rather 
low as pc ~ 5.5 X 10^^ g/cm^ and pc ^ 4.1 x 10^^ g/cm^, 
respectively. As we show in Sec. |IV[ the properties of the 
material ejected from the merger of binary neutron stars 
depend strongly on the radius of the neutron stars or pc. 



B. Initial conditions 

We employ binary neutron stars in quasiequilibria for 
the initial condition of numerical simulations as in our 
series of papers [25l [26]. The quasiequilibrium state is 
computed in the framework described in [39 to which 
the reader may refer. The computation of quasiequilib- 
rium states is performed using the spectral-method li- 
brary LORENE [41 . 

Numerical simulations were performed, systematically 
choosing wide ranges of the total mass and mass ratio of 
binary neutron stars. Because the mass of each neutron 
star in the observed binary systems is in a narrow range 
between ~ 1.2 - 1.45M0 we basically choose the 

neutron-star mass 1.20, 1.25, 1.30, 1.35, 1.40, 1.45, and 
1.5M0. Also, the mass ratio of the observed system q := 
mi/m2{< 1) where mi and m2 are lighter and heavier 
masses, respectively, is in a narrow range ~ 0.85 - 1. 
Thus, we choose as 0.8 < q < 1. Specifically, the 
simulations were performed for the initial data listed in 
Table [m 

The initial data were prepared so that the binary has 
about 3-4 quasicircular orbits before the onset of the 
merger. For four EOSs chosen, this requirement is ap- 
proximately satisfied with the choice of the initial an- 
gular velocity mQo = 0.026 for APR4 and ALF2 and 
ml^o = 0.025 for H4 and MSI. Here, m = mi + m2. 
For the following, the model is referred to as the name 
"EOS"-"mi" "m2"; e.g., the model employing APR4, 
mi = I.2M0, and m2 = 1.5M0 is referred to as model 



APR4-120150. 



III. FORMULATION AND NUMERICAL 
METHODS 

Numerical simulations were performed using an 
adaptive-mesh refinement (AMR) code SACRA [42] (see 
also [43] for the reliability of SACRA). The formulation, 
the gauge conditions, and the numerical scheme are ba- 
sically the same as those described in [42 , except for 
the improvement in the treatment of the hydrodynamics 
code for a far region. Thus, we here only briefiy review 
them and describe the present setup of the computational 
domain for the AMR algorithm and grid resolution. 



A. Formulation and numerical methods 

SACRA solves Einstein's evolution equations in the 
Baumgarte-Shapiro-Shibata-Nakamura formalism with a 
moving-puncture gauge [44 . It evolves a conformal fac- 
tor W := 7~^/^, the conformal three- metric jij := 
7~^^^7ij7 the trace of the extrinsic curvature a 
conformally- weighted trace-free part of the extrinsic cur- 
vature Aij := 'y~^^^{Kij—Kjij/3), and an auxiliary vari- 
able := —djY^. Introducing an additional auxiliary 
variable and a parameter 775, which we typically set 
to be ~ 0.8/m in units of c = G = Mq = 1, we employ a 
moving-puncture gauge in the form [45 

{dt-f3^dj)a = -2aK, (7) 
{dt-p^d,)p^ = {3/4)B\ (8) 
{dt-p'dj)B' = {dt-p'dj)r-risB\ (9) 

We evaluate the spatial derivative by a fourth-order cen- 
tral finite difference except for the advection terms, which 
are evaluated by a fourth-order lopsided upwind finite 
differencing scheme, and employ a fourth-order Runge- 
Kutta method for the time integration. 

To solve hydrodynamics equations, we evolve := 
pau^W~^^ Ui := hui^ and := hau^ — P/{pau^). Here, 
denotes the four velocity of the fiuid. The advection 
terms are handled with a high-resolution central scheme 
by Kurganov and Tadmor [46] with a third-order piece- 
wise parabolic interpolation for the cell reconstruction. 
We note that the total rest mass of the system is calcu- 
lated by 

= y p^d^x. (10) 

For the EOS employed in the numerical simulation, we 
decompose the pressure and specific internal energy into 
cold and thermal parts as 

P = ^cold + i^th , e = ^cold + £th' (11) 
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FIG. 2: Left: The gravitational mass as a function of the central density pc for spherical neutron stars in APR4, ALF2, H4, 
and MSI EOSs (the solid, dashed, dotted, and dash-dotted curves). Right: The same as the left panel but for the gravitational 
mass as a function of the circumferential radius. 



We calculate the cold parts of both var iables using the 
piecewise polytropic EOS (see section II A) from the 
primitive variable p, and then the thermal part of the spe- 
cific internal energy is defined from £ as £th = ^ — ^coid(p)- 
Because £th vanishes in the absence of shock heating, Sth 
is regarded as the finite temperature part determined by 
the shock heating in the present context. In this paper, 
we adopt a F-law ideal gas EOS for the thermal part as 



^th = (Fth - l)peth. 



(12) 



Following the conclusion of a detailed study in [47], Fth 
is chosen in the range 1.6 - 2.0 with the canonical value 
1.8. For several models, we performed simulations vary- 
ing the value of Fth , and explored the effects of the shock 
heating; as shown in Sec. |IV| numerical results depend 
fairly strongly on the value of Fth (although the depen- 
dence on Fth is not as strong as the dependence on the 
EOS, Pcoid). 

Because the vacuum is not allowed in any conservative 
hydrodynamics scheme (e.g., to derive the velocity by di- 
viding the momentum density by the density) , we put an 
artificial atmosphere of small density outside the neutron 
stars. The atmosphere has to be as tenuous as possible 
because a dense atmosphere may significantly affect the 
motion of the material ejected from binary neutron stars. 
Specifically, we set the density of the atmosphere in the 
following simple rule 



Pat 



J ^at Pmax ^ ^uni), 

I /at Pmax {r/r^ni) " {r> runi), 



(13) 



where Pmax denotes the maximum rest-mass density of 
the neutron stars at the initial state <, 10^^ g/cm^ (see 
Table |l|. We typically set /at = 10"^°, n = 3, and 
^uni = ISLmin whcrc 2Lmin dcuotcs the sidc length of 
the finest computational domain in the AMR algorithm 
(see Sec. IIIC and Table III). For MSI, a computational 
region is wider and we set /at = 10~^^ to reduce the at- 
mosphere mass. In these settings, the total rest mass of 



the atmosphere is always ~ 10 "^M© or less. In test sim- 
ulations, we also adopted n = 2 and /at = 10~^^ - 10~^^, 
and found that the numerical results on the ejected ma- 
terial such as its mass and its total energy depend only 
weakly on the values of n and /at (e.g., the ejected mass 
increases by ~ 10% if we change n from 3 to 2 (denser 
one) for some models of APR4 and H4). Hence, we could 
safely conclude that the tenuous atmosphere chosen in 
this work does not significantly affect the properties of 
the ejected material. 

We extracted I = \m\ = 2 modes of gravitational 
waves, hj^ and ^ hy calculating the outgoing part 
of the complex Weyl scalar ^4 at finite coordinate radii 
r ~ 2OOM0 - 4OOM0 and by integrating ^4 twice in time 
as in [33 , to which the reader may refer (see also [48]). 
We also analyzed the evolution of gravitational- wave fre- 
quency, which is determined by extracting the phase 
of ^4, arg(^4), and by taking the time derivative as 
27r/ := (i(arg(^4))/(it. To find the characteristic fre- 
quency of gravitational waves, we also define the average 
value of / by 



/ave • — 



j f\h\dt 
J \h\dt 



(14) 



where we used \h\ = {h'^ + ^x)"*^^^ weight factor. 

Then, we define the physical dispersion of / by 



{f-Uef\h\dt 



\h\dt 



(15) 



In the following, /ave and a/ are calculated for gravi- 
tational waves emitted by the remnant massive neutron 
stars. 
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TABLE II: List of the parameters of the initial condition for binaries chosen in numerical simulations: Total mass, mass ratio, 
masses of two components, initial value of angular velocity, and initial frequency of gravitational waves (/o = Qq/tt). 



Model 


m(M0) 


q 


mi (Mo) 


m2(M0 


) mQo 


/o (Hz) 


APR4-130160 


2.90 


0.813 


1.30 


1.60 


0.026 


579 


APR4-140150 


2.90 


0.933 


1.40 


1.50 


0.026 


579 


APR4-145145 


2.90 


LOOO 


1.45 


1.45 


0.026 


579 


APR4-130150 


2.80 


0.867 


1.30 


1.50 


0.026 


600 


APR4-140140 


2.80 


LOOO 


1.30 


1.50 


0.026 


600 


APR4-120150 


2.70 


0.800 


1.20 


1.50 


0.026 


622 


APR4-125145 


2.70 


0.862 


1.25 


1.45 


0.026 


622 


APR4-130140 


2.70 


0.929 


1.30 


1.40 


0.026 


622 


APR4-135135 


2.70 


LOOO 


1.35 


1.35 


0.026 


622 


APR4-120140 


2.60 


0.857 


1.20 


1.40 


0.026 


646 


APR4-125135 


2.60 


0.926 


1.25 


1.35 


0.026 


646 


APR4-130130 


2.60 


LOOO 


1.30 


1.30 


0.026 


646 


ALF2-140140 


2.80 


LOOO 


1.40 


1.40 


0.026 


600 


ALF2-120150 


2.70 


0.800 


1.20 


1.50 


0.026 


622 


ALF2-125145 


2.70 


0.862 


1.25 


1.25 


0.026 


622 


ALF2-130140 


2.70 


0.929 


1.30 


1.40 


0.026 


622 


ALF2-135135 


2.70 


1.000 


1.35 


1.35 


0.026 


622 


ALF2-130130 


2.60 


1.000 


1.30 


1.30 


0.026 


646 


H4-130150 


2.80 


0.867 


1.30 


1.50 


0.025 


577 


H4-140140 


2.80 


1.000 


1.40 


1.40 


0.025 


577 


H4-120150 


2.70 


0.800 


1.20 


1.50 


0.025 


598 


H4-125145 


2.70 


0.862 


1.25 


1.25 


0.025 


598 


H4-130140 


2.70 


0.929 


1.30 


1.40 


0.025 


598 


H4-135135 


2.70 


1.000 


1.35 


1.35 


0.025 


598 


H4-120140 


2.60 


1.000 


1.30 


1.30 


0.025 


621 


H4-125135 


2.60 


1.000 


1.30 


1.30 


0.025 


621 


H4-130130 


2.60 


1.000 


1.30 


1.30 


0.025 


621 


MSl-140140 


2.80 


LOOO 


1.40 


1.40 


0.025 


577 


MSl-120150 


2.70 


0.800 


1.20 


1.50 


0.025 


598 


MSl-125145 


2.70 


0.862 


1.25 


1.25 


0.025 


598 


MSl-130140 


2.70 


0.929 


1.30 


1.40 


0.025 


598 


MSl-135135 


2.70 


1.000 


1.35 


1.35 


0.025 


598 


MSl-130130 


2.60 


1.000 


1.30 


1.30 


0.025 


621 



B. Analysis of the ejected material 



In this section, we describe the method for analyzing 
the material ejected from the merger of binary neutron 
stars. Here, the ejected material is composed of a fluid 
element which is unbound by the gravitational potential 
of binary neutron stars and an object formed after the 
merger. Thus, first of all, we have to determine which 
fluid elements are unbound. To assess this point for all 
the fluid elements, we calculate u^t^ = Uf at each grid 
point. Here, is a timelike vector (1,0,0,0) which is a 
Killing vector at spatial infinity. If \ut\ > 1, we consider 
that the fluid element there is unbound. 



Then we calculate the total rest mass, total energy (ex- 
cluding gravitational potential energy), and total internal 



energy of the fluid element of |iAt| > 1 by 

M*esc = / p^d^x, (16) 

^tot,esc = / Tf_,^n^n!'^d^x 
J\ut\>i 

p^e^d^x, (17) 



/ 

J \u 



'\ut\>l 

Uesc = I p^ed^x, (18) 



/ , 



where T^^^ is the stress-energy tensor, 

T^y = phu^Uy + Pg^y, (19) 

and is the unit timelike hyper surface normal. We note 
that the total energy is not uniquely defined by £^tot,esc 
for dynamical spacetimes, and thus, the total energy de- 
fined here should be considered as an approximate mea- 
sure for it. We here choose this expression for simplicity. 
We then define kinetic energy approximately by 

^*esc -^tot,esc -^*esc ^esc- (^0) 
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TABLE III: The grid structure for the simulation in our AMR algorithm. Ax is the grid spacing in the finest-resolution domain 
with L being the location of the outer boundaries along each axis and Lmin = NAx. i?diam/Ax denotes the numbers of grid 
assigned inside the semi-major diameter of the lighter and heavier neutron stars in the finest level. Ao is the gravitational 
wavelength for the initial configuration. The last column shows the values of Fth employed. 



Model 


Ax (km) i^diair 




L (km) 


Lmin (km) 


Ao (km) 


Fth 


APR4-130160 


0.172 


102, 


96 


2636 


10.3 


518 


1.8 




APR4-140150 


0.167 


102, 


101 


2572 


10.0 


518 


1.8 




APR4-145145 


0.166 


102, 


102 


2550 


10.0 


518 


1.8 




APR4-130150 


0.172 


102, 


98 


2636 


10.3 


500 


1.8 




APR4-140140 


0.167 


102, 


102 


2572 


10.0 


500 


1.8 




APR4-120150 


0.172 


103, 


98 


2644 


10.3 


482 


1.6, 1.8 


, 2.0 


APR4-125145 


0.174 


102, 


100 


2665 


10.4 


482 


1.8 




APR4-130140 


0.170 


103, 


101 


2609 


10.2 


482 


1.8 




APR4-135135 


0.169 


102, 


102 


2601 


10.2 


482 


1.6, 1.8 


, 2.0 


APR4-120140 


0.174 


102, 


99 


2679 


10.5 


464 


1.8 




APR4-125135 


0.174 


102, 


100 


2665 


10.4 


464 


1.8 




APR4-130130 


0.171 


102, 


102 


2629 


10.3 


464 


1.8 




ALF2-140140 


0.195 


102, 


102 


3001 


11.7 


500 


1.8 




ALF2-120150 


0.200 


102, 


98 


3065 


12.0 


482 


1.8 




ALF2-125145 


0.199 


102, 


100 


3054 


11.9 


482 


1.8 




ALF2-130140 


0.198 


102, 


101 


3044 


11.9 


482 


1.8 




ALF2-135135 


0.195 


103, 


103 


3001 


11.7 


482 


1.8 




ALF2-130130 


0.199 


102, 


102 


3054 


11.9 


464 


1.8 




H4-130150 


0.222 


102, 


98 


3429 


13.4 


480 


1.8 




H4-140140 


0.219 


102, 


102 


3358 


13.1 


480 


1.8 




H4-120150 


0.228 


102, 


96 


3501 


13.7 


463 


1.6, 1.8 


, 2.0 


H4-125145 


0.226 


102, 


98 


3465 


13.5 


463 


1.8 




H4-130140 


0.223 


102, 


100 


3430 


13.4 


463 


1.8 




H4-135135 


221 


102, 


102 


3393 


13.3 


463 


1.6, 1.8 


, 2.0 


H4-120140 


0.230 


101, 


98 


3537 


13.8 


446 


1.8 




H4-125135 


0.227 


102, 


100 


3494 


13.6 


446 


1.8 




H4-130130 


0.223 


103, 


103 


3430 


13.4 


446 


1.8 




MSl-140140 


0.237 


103, 


103 


3644 


14.2 


480 


1.8 




MSl-120150 


0.249 


101, 


97 


3823 


14.9 


463 


1.8 




MSl-125145 


0.244 


102, 


99 


3751 


14.7 


463 


1.8 




MSl-130140 


0.244 


101, 


100 


3751 


14.7 


463 


1.8 




MSl-135135 


0.242 


102, 


102 


3715 


14.5 


463 


1.8 




MSl-130130 


0.244 


102, 


102 


3751 


14.7 


446 


1.8 





We found irrespective of models that T^esc is much (by 
about 1-2 orders of magnitude) larger than Uqsc- 

To approximately analyze the configuration of the 
ejected material, we also calculate the moments of in- 
ertia defined by 



Iii,esc = / p^{x''Yd?x, (no sum for i), (21) 
and then, define 



X 



Y 



(22) 



and R = x/X^ + F^. From dR/dt and dZ/dt, we can 
determine the typical (average) velocity of the ejected 
material, which is denoted by V^^^ and V^f^ in the follow- 
ing. 

We consider a model that the configuration of the 
ejected material is approximated by an axisymmetric 



anisotropic shell of uniform density as 

Pesc 7r/2-0o<e< 7r/2 + Oo 
p = ^ and R- <r< , 

otherwise, 



(23) 



where pesc, and are time- varying parameters. In 
this case, 

4:7r 

M,esc = Y/9esc(i?+ -i?-)sin6'o, (24) 



1R%-R^_ 
5Rl-R^_ 



R' 



(S-sin'^o), (25) 



Thus for an axial ratio, 

Z 



VR 



R' 



(26) 



(27) 



8 



sin 00 is calculated as 



1 



(28) 



Hence, from the axial ratio calculated for a numerical re- 
sult of the ejected material, we can approximately define 
the extent in the 6 direction; e.g., for r]R = 0.4 and 0.5, 
^0 ^ 40° and 51°, respectively. 

C. Setup of AMR grids 

An AMR algorithm implemented in SACRA can pre- 
pare a fine-resolution domain in the vicinity of compact 
objects as well as a sufficiently wide domain that covers 
a local wave zone. In the present study, we prepare ad- 
ditional domains wider than those used in our previous 
studies [26l |32l [33] , to follow the motion of the material 
ejected during the merger of binary neutron stars for a 
sufficiently long time (longer than 10 ms). 

The chosen AMR grids consist of a number of compu- 
tational domains, each of which has the uniform, vertex- 
centered Cartesian grids with (2A^ + 1, 2A/' + l, A^ + 1) grid 
points for (x, z) with the equatorial plane symmetry at 
z = 0. Since we chose that the grid spacing for three di- 
rections is identical, the shape of each AMR domain is a 
half cube. We chose = 60 for the best resolved runs 
in this work, and all the results shown in the following 
were obtained with this resolution. We also performed 
simulations with A = 40 and 50 (or 48) for several cho- 
sen models to check the convergence of the results (see 
Appendix A). 

We classify the domains of the AMR algorithm into 
two categories: one is a coarser domain, which covers a 
wide region including both neutron stars with its origin 
ffxed at the approximate center of mass throughout the 
simulation. The other is a ffner domain, two sets of which 
comove with two neutron stars and cover the region in 
their vicinity. We denote the side length of the largest 
domain, number of the coarser domains, and number of 
the ffner domains by 2L, and respectively. In 
this work, /c = 5 and = 4 (in total, 13). The grid 
spacing for each domain is hi = L/(2^A), where / = - 
^max(= /c + ^/ — 1) is the depth of each domain. In the 
following, we denote L/2^"^^^ by Lmin and hi^^^ by Ax. 

Table [lll| summarizes the parameters of the grid struc- 
ture for the simulations. For all the simulations, L is set 
to be L/c ^ 10 ms. This implies that the material cannot 
escape from the computational domain in ~ 10 ms after 
the onset of the merger, even if it could move with the 
speed of light. In reality, the speed of most of the ejected 
material is smaller than ~ 0.5c, and thus, the material 
stays in the second coarsest level for more than 10 ms. L 
is also much larger than the gravitational wavelengths at 
the initial instant Aq := tt/Qq. This implies that a spu- 
rious effect caused by outer boundaries when extracting 
gravitational waves is excluded in the present work more 
efficiently than in the previous works. The semi-major 



diameter of each neutron star is covered approximately 
by 100 grid points for A = 60. 



IV. NUMERICAL RESULTS 



Table IV summarizes the remnant formed, the rest 



mass and kinetic energy of the ejected material measured 
at 10 ms after the onset of the merger t = Emerge, and the 
characteristic (average) frequency of gravitational waves 
emitted by the hypermassive neutron star (HMNS) for 
A = 60 [61]. Here, t^erge is chosen to be the time at 
which the amount of the rest mass of the ejected mate- 
rial steeply increases. In the following two subsections, 
we summarize the results for the formation of HMNSs 
and black holes separately. 



A. Properties of the merger and mass ejection: 
HMNS case 



Binary neutron stars in quasicircular orbits evolve due 
to the gravitational-wave emission. Their orbital separa- 
tion decreases gradually, and eventually, the merger sets 
in. Previous studies (e.g., [26 ) clariffed that soon after 
the onset of the merger, either a long-lived HMNS or a 
black hole is formed. For most of the simulations in this 
paper performed with stiff EOSs and with the canoni- 
cal total mass 2.6 - 2.8M0, we found that a long-lived 
HMNS is formed with its lifetime much longer than its 
dynamical timescale ~ 0.1 ms and its rotation period 
~ 1 ms; the lifetime is longer than 10 ms for most of the 
models employed in this paper. In this section, we pay 
attention to the case that such a HMNS is formed. 

Figures [3] - [5] display snapshots of the density proffles 
in the merger for models APR4-135135, APR4-120150, 
and H4-120150, respectively. Figure [6| also displays the 
central density as a function of time for the models with 
mi = m2 = I.35M0 (left), and mi = I.2M0 and 
m2 = I.5M0 (right). These ffgures show that a compact 
and nonaxisymmetric object (proto HMNS) is formed in 
the central region soon after the onset of the merger. The 
shape and compactness of the HMNS depend strongly on 
the EOS and mass ratio; e.g., the presence of the asym- 
metric spiral arms found in the top panels of Figs. |4] and [5] 
is the feature only for the asymmetric binaries; the ampli- 
tude of the quasiradial oscillation is larger for the equal- 
mass binaries; a high-amplitude quasiradial oscillation 
is a unique property found only for models with APR4 
(see Fig. [6]). However, it is universal that the HMNSs are 
rapidly rotating and nonaxisymmetric, irrespective of the 
EOS, total mass (m < 2.8M0), and mass ratio, as found 
in previous studies [211 ISSl [26] . This rapid rotation to- 
gether with the nonaxisymmetric conffguration not only 
results in the emission of strong gravitational waves but 
also is the key for an efficient mechanism of angular mo- 
mentum transport from the HMNS to the surrounding 
material because the HMNS exerts the torque. 
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FIG. 3: Snapshots of the density profile for the merger of binary neutron stars for an equal-mass model APR4-135135. The 
first row shows the density profiles in the equatorial plane and in the central region, and second - fourth ones show the density 
profile for a wide region in the x-y, x-z, and y-z planes, tmerge ~ 11.3 ms for this model. 
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FIG. 4: The same as Fig. |3] but for unequal-mass model APR4-120150. tmerge ~ 10.3 ms for this model. 
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FIG. 5: The same as Fig.[4]but for models H4-120150. tmerge ~ 8.8 ms for this model. 
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FIG. 6: The central density as a function of time for models with mi — m2 — 1.35M0 (left), and mi = 1.2M0 and m2 = 1.5Mq 
(right). Before the merger of unequal mass binaries, the central density of heavier neutron stars are plotted. Fth = 1.8 is 
employed for the results presented here. 
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FIG. 7: Snapshots of the thermal part of the specific internal energy (sth) profile in the vicinity of HMNSs on the equatorial 
(top) and x-z (bottom) planes for an equal-mass model APR4-135135. The rest-mass density contours are overplotted for every 
decade from 10^^ g/cm^. 



Figures [3]-|5] indicate that there are two important pro- 
cesses for the mass ejection. The first one is the heating 
by shocks formed at the onset of the merger between the 
inner surfaces of two neutron stars. Figures [7| and [8] dis- 
play snapshots of the thermal part of the specific internal 
energy, Sth, in the vicinity of HMNSs for APR4-135135 



and APR4-120150, respectively. These figures show 
clearly that hot materials with ^th ^0.1 (^ lOOMeV) are 
indeed ejected from the HMNSs, in particular, to bidirec- 
tional regions on the equatorial plane and to the polar 
region. This suggests that the shock heating works effi- 
ciently to eject materials from the HMNSs. This occurs 
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FIG. 8: The same as Fig. [S] but for an unequal-mass model APR4-120150. 



in an outstanding manner in particular for the equal-mass 
(and only slightly asymmetric) binaries. The heated- up 
material is pushed outwards by the thermal pressure gen- 
erated by the shock approximately in the plane parallel 
to the (rotating) shock surface. Subsequently, it expands 
outwards with rotation, and eventually forms hot spiral 
arms around the HMNS. This component subsequently 
gains angular momentum (and hence kinetic energy) due 
to the torque exerted by the HMNS of a nonaxisymmetric 
configuration, and a fraction of the material eventually 
gains the kinetic energy that is large enough for it to es- 
cape from the system. This effect plays a primary role 
for the early mass ejection that occurs in the first a few 
ms after the onset of the merger. 

A stronger shock appears to play basically a positive 
role for increasing the amount of the ejected material, be- 
cause the amount of the heated- up material can be more, 
and as a result, the materials in the spiral arm and ejected 
fraction increase. A stronger shock is formed for softer 
EOSs or for binaries composed of more compact neutron 
stars (e.g., APR4 in the context of canonical-mass neu- 
tron stars). The reason is that neutron stars for such an 
EOS can achieve a more compact state (cf. Fig. [6| and 
at the merger, the collision velocity of two neutron stars 
is larger (the minimum separation between two stars is 
smaller). This point will be in more detail described in 
Sec. |IV A 3 A strong shock could be also formed for bi- 



naries with the total mass close to the critical value for 
the collapse to a black hole even for stiff EOSs, because a 
highly compressed state is realized by the strong gravity. 

The shocks are also formed continuously in the outer 
part of the HMNS during its evolution through the in- 
teraction with spiral arms formed in its envelope due to 
a torque exerted by the HMNS (see below). This effect 
plays an important role in a relatively longer-term mass 
ejection with the duration ~ 10 - 20 ms. 

The secondly important process for the mass ejection 
is a hydrodynamic interaction induced by the HMNS of 
a nonaxisymmetric configuration that exerts the torque 
to the surrounding material and transports the angular 
momentum outwards. Since it is rapidly rotating, the 
HMNS works as an efficient torque supplier. Our simula- 
tions show that this process is important in particular in 
the early phase of the merger: For the nearly equal-mass 
binaries, a fraction of the material that spreads outwards 
by the shock formed at the merger subsequently gains an- 
gular momentum from the HMNS and eventually obtains 
kinetic energy large enough to escape from the system; 
for sufficiently asymmetric binaries (for small values of 
g), a less-massive neutron star is tidally elongated dur- 
ing the early phase of the merger, a fraction of its mate- 
rial forms spiral arms, and it subsequently gains angular 
momentum from the HMNS enough to escape from the 
system. In the early mass ejection caused by the torque 
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exerted by the HMNS, the material is primarily ejected 
in the direction near the equatorial plane, and the typi- 
cal velocity of the escaping material in this early stage is 
quite high ~ 0.5 - 0.8c (follow the locations of the head 
of the ejected materials in Figs. [3]-|5|. The maximum 
velocity is larger for the EOS that yields smaller-radius 
neutron stars; for APR4, it is ~ 0.8c and for MSI, it is 
~ 0.5c. This also depends on the mass ratio for models 
with a large neutron-star radius (for models of H4 and 
MSI). 

In the later phase, the mass ejection appears to oc- 
cur by the combination of the shock heating and by the 
torque exerted by the HMNS. As mentioned already, the 
continuous shock heating occurs in the envelope of the 
HMNS in the presence of spiral arms. Due to this, a 
fraction of the material gains large kinetic energy. In ad- 
dition, the material in the outer region gains angular mo- 
mentum by the torque exerted by the HMNS. These two 
effects give a fraction of the material the escape velocity. 
By this process, the material is gradually ejected from 
the system in a quasispherical manner; the anisotropy of 
the configuration of the ejected material is not as large as 
that of the material ejected in the early stage. This indi- 
cates that the shock heating plays a relatively important 
role. The average velocity of the escaping material in this 
process is sub-relativistic ~ 0.15 - 0.25c (see Table IV). 



In the mass ejection process, these two nonlinear ly cou- 
pled effects (shock heating and torque exerted by the 
HMNS) play a substantial role. As a result, the amount 
of the ejected material depends on the EOS, the total 
mass of the system, and the mass ratio in a nonlinear 
manner. Thus, a small change (associated, e.g., with the 
grid resolution, the initial orbital separation, configura- 
tion of the atmosphere, and presence or absence of the tt 
symmetry for equal-mass binaries) results in the change 
in the rest mass and kinetic energy of the ejected ma- 
terial; this fluctuation is in general small, ~ 10 - 20%, 
for unequal-mass binaries for which the torque plays a 
primary role (see Appendix A). For the equal-mass case, 
the convergence is poor because a strong shock often oc- 
curs at the merger and plays a primary role in the mass 
ejection. The possible reason for this poor convergence 
is that shocks are always computed by the first-order ac- 
curacy in the spatial grid resolution, and hence, the ac- 
curacy is low and in addition, the ejected mass is a tiny 
part of the entire system. A random error for the entire 
system computed with a low accuracy significantly (and 
randomly) affects a tiny part (i.e., the ejected material), 
resulting in the poor convergence. (We note that for 
global quantities, the convergence is usually good.) For 
some models (such as ALF2-135135 and MSl-135135), 
the ejected mass increases steeply with the grid resolu- 
tion, and for such cases, the results in this paper might 
give the lower bound. 

In the following subsections, we describe the properties 
of the ejected material in more detail. 



1. Gravitational waves 

First of all, we summarize the properties of grav- 
itational waves emitted by the HMNS, because its 
gravitational- wave frequency, which is determined by the 
spin of the HMNS, has a correlation with the amount of 
the ejected material 

As mentioned already, HMNSs exert the torque to its 
surrounding material. The efficiency of the angular mo- 
mentum transport is higher, in general, for the faster 
rotating and more compact HMNS. Associated with this 
property, the frequency of gravitational waves and the ef- 
ficiency of the angular momentum transport are expected 
to be closely related. The characteristic spin frequency 
for these deformed HMNSs can be determined from grav- 
itational waves emitted by them. Figure [9] displays grav- 
itational waves and their frequency as functions of time 
for eight models with mass (mi, 7712) = (1.3M0, 1.4M0) 
and (mi,m2) = (I.2M0, 1.5M0) and with four EOSs. 
These plots show that quasiperiodic gravitational waves 
are emitted by the HMNSs for all the models. Namely, 
the gravitational- wave frequency does not change signifi- 
cantly during the evolution of the HMNSs. However, the 
frequency is not constant exactly and actually varies with 
time. This is natural because (i) the HMNSs quasiradi- 
ally oscillate with time in their early stage of the evolu- 
tion, and (ii) the HMNSs lose the energy and angular mo- 
mentum due to the gravitational- wave emission and hy- 
drodynamic angular momentum transport process, and 
hence, their configuration evolves. These two effects re- 
sult in the variation in the characteristic spin velocity 
and frequency of gravitational waves. The degree of the 
variation in the frequency of gravitational waves is larger 
(a) for the HMNS with the EOS that yields a compact 
neutron star (we often call such an EOS soft EOS in 
this paper), and (b) for the HMNS for which the mass 
is close to the critical value to the collapse to a black 
hole; see, e.g., the gravitational- wave frequency for model 
ALF2-130140. The case (a) is due to the fact that at the 
merger, the central density significantly increases in the 
soft EOSs, resulting in subsequent high- amplitude oscil- 
lations. The case (b) is due to the fact that for such a 
HMNS, a small change in the spin velocity results in a 
large change in the central density. 

Figure [lO] plots the Fourier spectra for gravitational 
waves shown in Fig. |9] This shows that there are peaks 
for a high-frequency band 2 kHz ^ / ^ 4 kHz irrespec- 
tive of models. For a "soft" EOS that yields a compact 
neutron star for the canonical mass, the peak frequency 
is higher (e.g., for the spectra of APR4, the peak fre- 
quency is the highest among the four EOSs), and a cer- 
tain correlation exists between the peak frequency and 
stellar radius [49 . The peak frequency is approximately 
associated with the typical frequency of quasiperiodic os- 
cillation of gravitational waves found in Fig. [9j How- 
ever, as already mentioned, the (nonaxisymmetric) oscil- 
lation frequencies of the HMNSs vary during the evolu- 
tion due to a quasiradial oscillation and the back reaction 
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FIG. 9: Gravitational waves (h-^D/m) and the frequency of gravitational waves / as functions of time for models APR4-130140 
(top left), APR4-120150 (top right), ALF2-130140 (second top left), ALF2-120150 (second top right), H4-130140 (third left), 
H4-120150 (third right), MSl-130140 (bottom left), and MSl-120150 (bottom right). For ALF2-130140, a black hole is formed 
at 11 ms after the onset of the merger, and ringdown gravitational waves are emitted in the final phase. For all the panels, the 
vertical axis shows the non-dimensional amplitude, /i+D/m, with D being the distance to the source. 
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TABLE IV: Summary of numerical results. The remnant, the total rest mass, M*esc, the kinetic energy, T*esc, the R and Z 
components of the average velocity of escaping material, and t^fc, of the ejected material, and characteristic frequencies of 
gravitational waves emitted by HMNSs for 5 and 10 ms time integration after the formation of the HMNSs. The total rest mass, 
kinetic energy, and average velocity are measured at ^ 10 ms after the onset of the merger. The dispersion of /ave shown here is 
CT/. BH denotes black hole. The remnant is judged at ~ 30 ms after the onset of the merger. All the results shown are those in 
the run with = 60 and our standard setting of atmosphere. The rest mass and kinetic energy of the ejected material have the 
uncertainty of order 10%. The approximate lifetime of HMNSs for APR4-130150, APR4-140140, ALF2-140140, ALF2-130140, 
ALF2-135135, H4-130150, H4-140140, H4-135135 (Fth = 1.6)), and H4-135135 (Fth = 1.8)) is - 30, 30, 5, 10, 15, 20, 10, 15, 
25 ms for A" = 60, respectively. 



Model 


Fth 


Remnant 


M*esc(lO-^M0) 


T*esc(10^°ergs) 


vJ^c/c 




/ave, 5ms (kHz) 


/ave,10ms (kHz) 


APR4-130160 


1.8 


BH 


2.0 


1.5 


0.24 


0.08 


— 


— 


APR4-140150 


1.8 


BH 


0.6 


0.9 


0.35 


0.12 








APR4-145145 


1.8 


BH 


0.1 


< 0.1 


0.29 


0.13 








APR4-130150 


1.8 


HMNS^BH 


12 


8.5 


0.23 


0.12 


3.48 ± 0.47 


3.46 ± 0.37 


APR4-140140 


1.8 


HMNS^BH 


14 


10 


0.22 


0.15 


3.53 ±0.52 


3.52 ±0.48 


APR4-120150 


1.6 


HMNS 


9 


5 


0.20 


0.10 


3.47 ±0.30 


3.44 ±0.27 


APR4-120150 


1.8 


HMNS 


8 


5.5 


0.23 


0.11 


3.44 ±0.30 


3.41 ±0.24 


APR4-120150 


2.0 


HMNS 


7.5 


5.5 


0.24 


0.12 


3.32 ±0.32 


3.27 ±0.26 


APR4-125145 


1.8 


HMNS 


7 


4.5 


0.22 


0.11 


3.36 ±0.31 


3.31 ±0.25 


APR4-130140 


1.8 


HMNS 


8 


5 


0.19 


0.12 


3.30 ±0.29 


3.27 ±0.28 


APR4-135135 


1.6 


HMNS 


11 


6 


0.19 


0.13 


3.46 ± 0.42 


3.45 ±0.37 


APR4-135135 


1.8 


HMNS 


7 


4 


0.19 


0.12 


3.31 ±0.35 


3.31 ±0.32 


APR4-135135 


2.0 


HMNS 


5 


3 


0.19 


0.13 


3.35 ±0.39 


3.33 ±0.33 


APR4-120140 


1.8 


HMNS 


3 


2 


0.21 


0.12 


3.15 ±0.21 


3.13 ±0.19 


APR4-125135 


1.8 


HMNS 


5 


3 


0.18 


0.10 


3.22 ±0.25 


3.19 ±0.24 


APR4-130130 


1.8 


HMNS 


2 


1 


0.19 


0.10 


3.22 ±0.28 


3.19 ±0.26 


ALF2-140140 


1.8 


HMNS^BH 


2.5 


1.5 


0.21 


0.13 


2.93 ±0.42 


— 


ALF2-120150 


1.8 


HMNS 


5.5 


3 


0.21 


0.10 


2.70 ±0.19 


2.71 ±0.16 


ALF2-125145 


1.8 


HMNS 


3 


1.5 


0.20 


0.10 


2.66 ±0.14 


2.66 ±0.13 


ALF2-130140 


1.8 


HMNS ^ BH 


1.5 


0.8 


0.16 


0.11 


2.73 ±0.19 


2.75 ±0.17 


ALF2-135135 


1.8 


HMNS ^ BH 


2.5 


1.5 


0.22 


0.12 


2.75 ±0.18 


2.76 ±0.16 


ALF2-130130 


1.8 


HMNS 


2 


1.0 


0.19 


0.10 


2.58 ±0.18 


2.56 ±0.16 


H4-130150 


1.8 


HMNS^BH 


3 


2 


0.19 


0.10 


2.44 ±0.17 


2.45 ±0.15 


H4-140140 


1.8 


HMNS^BH 


0.3 


0.2 


0.17 


0.13 


2.63 ±0.23 


2.77 ±0.41 


H4-120150 


1.6 


HMNS 


4.5 


2 


0.19 


0.10 


2.28 ±0.16 


2.29 ±0.14 


H4-120150 


1.8 


HMNS 


3.5 


2 


0.21 


0.09 


2.30 ±0.18 


2.31 ±0.15 


H4-120150 


2.0 


HMNS 


4 


2 


0.21 


0.09 


2.24 ±0.15 


2.23 ±0.14 


H4-125145 


1.8 


HMNS 


2 


1.5 


0.19 


0.10 


2.41 ±0.15 


2.41 ±0.13 


H4-130140 


1.8 


HMNS 


0.7 


0.4 


0.18 


0.10 


2.42 ±0.17 


2.42 ±0.15 


H4-135135 


1.6 


HMNS^BH 


0.7 


0.4 


0.21 


0.11 


2.49 ±0.19 


2.54 ±0.16 


H4-135135 


1.8 


HMNS^BH 


0.5 


0.2 


0.19 


0.11 


2.44 ±0.20 


2.48 ±0.16 


H4-135135 


2.0 


HMNS 


0.4 


0.2 


0.20 


0.10 


2.39 ±0.21 


2.43 ±0.17 


H4-120140 


1.8 


HMNS 


2.5 


1 


0.19 


0.10 


2.30 ±0.15 


2.30 ±0.14 


H4-125135 


1.8 


HMNS 


0.6 


0.3 


0.18 


0.10 


2.29 ±0.17 


2.27 ±0.14 


H4-130130 


1.8 


HMNS 


0.3 


0.1 


0.16 


0.10 


2.35 ±0.18 


2.38 ±0.14 


MSl-140140 


1.8 


MNS 


0.6 


0.2 


0.13 


0.09 


2.09 ±0.14 


2.06 ±0.12 


MSl-120150 


1.8 


MNS 


3.5 


1.5 


0.19 


0.10 


2.08 ±0.11 


2.09 ±0.09 


MSl-125145 


1.8 


MNS 


1.5 


0.8 


0.19 


0.11 


2.02 ±0.14 


1.99 ±0.15 


MSl-130140 


1.8 


MNS 


0.6 


0.2 


0.17 


0.09 


2.05 ±0.14 


2.02 ±0.13 


MSl-135135 


1.8 


MNS 


1.5 


0.6 


0.14 


0.08 


1.98 ±0.18 


1.95 ±0.16 


MSl-130130 


1.8 


MNS 


1.5 


0.5 


0.15 


0.08 


1.93 ±0.19 


1.90 ±0.17 



due to the gravitational-wave emission and angular mo- 
mentum transport process, and hence, the peak frequen- 
cies change with time, resulting in the broadening of the 
peak or appearance of the multi peaks. Therefore, it is 
not a very good idea to determine the characteristic fre- 
quency from the peak of the Fourier spectrum. Rather, 
the Fourier spectrum might provide an inaccurate mes- 
sage when we determine the characteristic oscillation fre- 



quency. Thus, we determine the average frequency from 
the results of the frequency shown in Fig. |9] in terms 
of Eq. (14) with the dispersion determined by Eq. (15). 



Here, the time integration is performed for 5 and 10 ms 
after the formation of the HMNSs. The last two columns 



of Table [IV| list the average frequency and the disper- 
sion determined for 5 and 10 ms integration. Note that 
the typical nonaxisymmetric oscillation frequency of the 
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FIG. 10: Fourier spectra of gravitational waves for the results shown in Fig. [5] The amplitude is shown for the hypothetical 
event at a distance of 50 Mpc along the direction perpendicular to the orbital plane (the most optimistic direction) . The black 
dot-dot curve is the noise spectrum of the advanced LIGO with an optimistic configuration for the detection of high-frequency 
gravitational waves (see https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=2974). 



HMNSs is half as large as the values listed in Table IV 
because the listed ones are the gravitational-wave fre- 
quencies. 

The value of the oscillation frequency for a given mass 
of the HMNS depends primarily on its radius, i.e., a 
stiffness of the EOS. For the EOS that yields small- 
radius neutron stars ("soft" EOS), the oscillation fre- 
quency and peak frequency of gravitational waves are 
higher, because the spin angular velocity of the HMNS 
is close to the Kepler velocity, and thus, the oscillation 
and peak frequencies are qualitatively proportional to 
(Mhmns/^hmns)^^^ where Mhmns and i^HMNS denote 
the typical mass and radius of a HMNS. The oscillation 
frequency depends also weakly on the value of Tth^ For 
the smaller value of it, the frequency is slightly higher for 
many cases, because the effect of shock heating is weaker, 
and the HMNS becomes more compact. 

For a larger spin of the nonaxisymmetric HMNS, the 
material surrounding the HMNS can receive a torque 
with a higher efficiency. This suggests that for the merger 
of a binary neutron star composed of smaller-radius neu- 
tron stars, the amount of the ejected material could be 
larger. As shown in Sec. IV A 3[ this is indeed the case 
(in particular for unequal-mass models), as long as the 
models in this paper are concerned. 

Table |IV| as well as Fig. [lO] also show that the magni- 
tude of the dispersion, crj, is not negligible. For APR4 
for which the neutron-star radius is rather small and the 
amplitude of a quasiradial oscillation induced at the for- 
mation of the HMNSs is rather large, the magnitude of 
the dispersion is 0.2 - 0.5 kHz; for m = 2.7M0, the typi- 
cal value is 0.3 - 0.4 kHz. For other EOSs, the dispersion 
is relatively small. However, it is still 0.1-0.2 kHz. Thus, 
we conclude that the characteristic frequency of gravita- 
tional waves from HMNSs varies with time in general. 
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FIG. 11: R{= (i?^)^/^) and Z(= {z^)^^^) as functions of 
time for APR4-120150, ALF2-120150, H4-120150, and MSl- 
120150. 



2. Average velocity of the ejected material 

Figure pT] plots R and Z as functions of time for APR4- 
120150, ALF2-120150, H4-120150, and MSl-120150. 
Note that dR/dt and dZ/dt may be considered as an 
average velocity of the ejected material in the cylin- 
drical and vertical directions, respectively, and that a 
similar result is found for other choices of mass. This 
shows that the material ejected expands with an approx- 
imately constant sub-relativistic velocity ~ 0.15 - 0.25c 
for t — ^merger ^ 2 ms in the cylindrical direction and 
the velocity in the vertical direction is 0.4 - 0.5 times as 
large as that in the cylindrical direction. This suggests 
that the vertical thickness angle of the ejected material, 
6>o, is ~ 40 - 50°. Namely, the ejected material expands 
in a moderately anisotropic manner. Note that the veloc- 
ity in the cylindrical direction is primarily caused by the 
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torque exerted by the HMNS, while the velocity in the 
vertical direction is primarily caused by the shock heat- 
ing. This implies that both effects play an important 
role. 

The velocity in the later phase, t — tmerge ^ 3 ms, is 
larger for APR4 than that for other EOSs employed in 
this paper. This is due to the fact that with APR4, a 
more compact state is realized in the HMNS, and hence, 
(i) a strong shock associated with the compression by 
a strong gravity and a subsequent large- amplitude oscil- 
lation (cf. Fig. [6| occurs, resulting in an efficient mass 
ejection, and (ii) the HMNS strongly exerts the torque to 
its surrounding material. For APR4, a relatively dense 
atmosphere surrounding the HMNS is formed not only in 
the vicinity of the equatorial plane but also in the vertical 
direction (compare Figs. [4] and|5|. This also reflects the 
fact that a strong shock heating occurs with this EOS 
(see Figs. [7| and 



Dependence on EOS 



Figure [12] plots the total rest mass and kinetic energy 
of the material ejected from the HMNSs as functions of 
t — tmerge for scvcral modcls; for the left and right pan- 
els, the masses of two neutron stars are (l.SM©, I.4M0) 
and (I.2M0, 1.5M0), respectively, with the total mass 
2.7M0, while four EOSs are chosen. This shows that the 
rest mass and kinetic energy of the ejected material de- 
pend strongly on the EOS. The primary reason is that the 
compactness of the HMNS depends strongly on the EOS. 
For APR4 and ALF2, neutron stars of canonical masses 
1.2 - I.5M0 have a relatively small radius (cf. Table 
This implies that the merger sets in at a compact orbit, 
and the formed HMNS is more compact than that formed 
in stiffer EOSs that yield large-radius neutron stars. 

A high compactness of a HMNS affects the properties 
of the material ejected from it in the following two ways. 
First, the HMNS is more rapidly rotating, and hence, it 
exerts the torque, caused by its nonaxisymmetric config- 
uration and rapid rotation, to the material in the outer 
region more efficiently than a less compact HMNS. As a 
result of this effect, a fraction of the material that gains 
the kinetic energy large enough to escape from the sys- 
tem is increased. In addition, during the formation of 
such a compact HMNS, a quasiradial oscillation with a 
high amplitude is often induced (see Fig. [6|. This is 
in particular the case for APR4 in which the EOS be- 
comes stiff for a high-density region although it is rather 
soft for the density of canonical-mass neutron stars. This 
quasiradial oscillation helps the material surrounding the 
HMNS to obtain k inetic energy through shock heating 
(see section IV A 1 ). 

A possibly important fact to be noted is that the mate- 
rial, which eventually escapes from the system, initially 
stays in the vicinity of the HMNS. Namely, this mate- 
rial stays in a deep gravitational potential well initially, 
and thus, it is trapped. For a more compact HMNS, 



this potential should be deeper, and hence, the material 
there needs to obtain more energy to escape from the 
HMNS. At the same time, however, such a material can 
gain a stronger torque and thermal energy for a longer 
timescale, because it is trapped for a longer duration, and 
as a result, the material could get more kinetic energy if 
the HMNS is more compact. 

Namely, there are two competing effects, and it is not 
trivial at all which effects are more important. If the 
trapping effect due to the deep potential well plays a 
more important role, we should find the evidence that 
less material is ejected from more compact HMNS. How- 



ever, Fig. 12 shows that the ejected rest mass is smaller 
for EOS with larger neutron star radii (less compact neu- 
tron star). This indicates that the trapping mechanism 
is less important than the effects of the quasiradial oscil- 
lation and the torque exerted from the HMNS, as long 
as the comp arison a mong four model EOSs is concerned 
(but see Sec. IV A 5 for an evidence that this may not be 
always the case). 

The compactness of HMNSs, and hence, the EOS of 
neutron stars, are well reflected in the frequency of grav- 
itational waves emitted by the HMNS, as already de- 



scribed in Sec. |IV A 1[ For a given total mass and mass 
ratio of the binary system, the frequencies are higher 
for binaries composed of more compact neutron stars 
("softer" EOS), because the formed HMNS is more com- 
pact and hence the rotational angular velocity approx- 
imately proportional to (Mhmns/^hmns)^^^ larger. 
Since we found that the rest mass and kinetic energy of 
the ejected material are larger for the EOS that yields 
more compact HMNSs, these quantities and the fre- 
quency of gravitational waves should have a correlation. 



Figure [13] plots the rest mass of the ejected material as 
a function of the characteristic frequency of gravitational 
waves emitted by the HMNS, /ave^ foi" several models. 
Here, /ave is determined by the 5 ms integration using 
Eq. ( [Ts] ). For this plot, the results with Fth = 1-8 are 
adopted. Note that the typical frequency for this plot 
is determined primarily by the chosen EOS. This figure 
shows that for a given mass ratio these two quanti- 
ties have a correlation; the total rest mass of the ejected 
material increases with the gravitational- wave frequency. 

Howev er, as al ready noted, there are counter examples 
(see Sec. IV A 5). Namely, for some cases, the shallow 



potential helps in enhancing the mass ejection. For such 
models, the correlation like that found in Fig. 13 does 
not hold. 



4. Dependence on Fth 

The total rest mass and kinetic energy for the ejected 
material depend also on the value of Fth- The possible 
reason is described as follows. 

For larger values of Fth, the effect of shock heating 
is stronger. This implies that the thermal energy of the 
HMNS is increased via the shock heating more efficiently. 



19 



0.01 
0.001 



S 0.0001 



10 



0.1 



APR4-130140 - 

ALF2-130140 - 

H4-130140 ■■ 
MSl-130140 

H > > > > \- 







t-Wrge (ins) 



10 



^ 0.01 

^ 0.001 

J 0.0001 
10 



0.1 



APR4-120150 
ALF2-120150 
H4-120150 
MSl-120150 



M. X:: 



10 



t - Wrge (ins) 



FIG. 12: M*esc and T*esc as functions of t - tmerge (left) for models APR4-130140, ALF2-130140, H4-130140, MSl-130140, and 
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gravitational- wave frequency emitted by the HMNS. Here, the 
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and thus, the material located outside the HMNS that 
will eventually escape from the system expands more ef- 
ficiently at the merger and during the subsequent shock 
heating. This effect could result in increasing the ejected 
material. 

On the other hand, the HMNS becomes less compact 
by more efficient shock heating for the larger value of Fth, 
and hence, the amplitude of the quasiradial oscillation is 
smaller. This suggests that although the outward veloc- 
ity of the material caused by the shock heating is initially 
larger for the larger values of Fth , the subsequent gain of 
the kinetic energy via the shock heating could be smaller. 
The less compact HMNS could be also less favorable for 
exerting the torque to its surrounding material because 
the rotational velocity is slower. Therefore, the total rest 
mass and kinetic energy of the material ejected from the 
system depend on two competing nonlinear processes, as 
in a mechanism similar to that mentioned in Sec. IIV A 31 



Figure 14 compares the evolution of R and Z for 
Fth = 1.6, 1.8, and 2.0 for models _APR4-135135 and 
APR4-120150. For APR4-135135, R is larger for the 
larger values of Fth for t — t 



merge 



^ 1.5 ms. This 
agrees with the prediction that the shock heating effect 
is stronger and the material expands in a wider region for 
the larger values of Fth- For Z, the similar result is found 
for 3 ^ t — tmerge ^ 5 ms. Howcvcr, for t — tmerge ^ 5 ms, 
dR/dt and dZ/dt have a similar magnitude depending 
only weakly on the value of Fth- This is due to the fact 
that the mass ejection is primarily driven by the torque 
exerted by the HMNS. 

For APR4-120150, soon after the onset of the merger, 
R and Z are only slightly larger for the larger values of 
Fth- This is due to the fact that the mass ejection is 
primarily driven by the tidal effect caused by the mass 
asymmetry irrespective of the values of Fth- However, for 
t — tmerge ^ 3 ms, dR/dt and dZ / dt become smaller for 
the smaller values of Fth- This is due to the fact that for 
the smaller value of Fth, in particular, for Fth = 1-6, new 
materials with a smaller velocity are gradually ejected 
for the later time. This occurs due to the fact that for 
the lower value of Fth, a longer-term mass ejection driven 
also by the torque exerted by the HMNS occurs. Namely, 
for both equal-mass and unequal-mass cases, a longterm 
mass ejection driven by the angular momentum transport 
from the HMNSs play an important role for Fth = 1-6. 



Figure [15] compares the evolution of the total rest mass 
and kinetic energy for Fth = 1.6, 1.8, and 2.0 for mod- 
els APR4-135135, APR4-120150, H4-135135, and H4- 
120150. All the panels of Fig. [15] clearly show that for 



the early time, t — tmerge ^ 1-5 ms for APR4-135135 and 
H4-135135, ^ 5 ms for APR4-120150 and H4-120150, 
these two quantities are larger for the larger value of Fth- 
Namely the stronger shock heating associated with the 
larger value of Fth plays an important role. However, af- 
ter the early time, the rest mass tends to be larger for 
the smaller value of Fth- In particular, for Fth = 1-6, a 
rapid increase in the total rest mass is found. Thus, a 
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CD 

o 
in 
O 



0.012 
0.01 
0.008 
^ 0.006 
I 0.004 
2 0.002 

6 
5 
4 
3 
2 
1 




^ 0.001 
0.0008 
0.0006 
0.0004 

g"" 0.0002 


^ 0.4 

5-1 

o 0.3 
o 

C 0.2 

I 0.1 
H 






t-Wrge (ins) 




^ 0-01 
0.008 
0.006 
0.004 
0.002 

if 
10 



O 
in 
O 



- -f/ 








^ ' ' ^ 


— 1 1 1 1 1 1 1 — ; 

r,h=2.o 1 


r ij ;?\ 


r,h=i-8 \ 




r,h=i-6 ; 











5 10 



t - Wrge (ins) 




t-Wrge (ins) 



t - Wrge (ms) 



FIG. 15: M*esc and T*esc as functions of t-Uerge (left) for models APR4-135135 (top left), APR4-120150 (top right), H4-135135 
(bottom left), and H4-120150 (bottom right) with Fth = 2.0, 1.8, and 1.6. 



longterm mass ejection process driven by the torque ex- 
erted by the HMNSs works for the smaller values of Fth 
(i.e., for more compact HMNSs), and this mechanism is 
remarkable for Fth = 1.6. 

For APR4-120150 and H4-120150, the rest mass of the 
ejected material is largest for Fth = 1.6. However, the ki- 
netic energy depends weakly on the value of Fth. This im- 



plies that although more materials are ejected, the gained 
kinetic energy is not very large for Fth = 1.6, because the 
velocity of material ejected later by the tidal torque is not 
very large. 

The dependence of the rest mass and kinetic energy 
of the ejected material on Fth is qualitatively similar for 
APR4 and H4. This indicates that the properties sum- 
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marized in this subsection would hold irrespective of the 
EOS. 



5. Dependence of the ejected material on the mass ratio 
and total mass 



The total rest mass and kinetic energy of the material 
ejected from the HMNSs depend also on the mass ratio 
and total mass of binary neutron stars. The degree of the 



dependence depends on the EOS. Figure [16] plots M^esc 
and T^esc as functions of t — Emerge for APR4 and H4 
with three mass ratios and with the total mass 2.1 Mq. 
For the models with H4, the total rest mass and kinetic 
energy of the ejected material depend strongly on the 
mass ratio; e.g., the total rest mass and kinetic energy 
for q = 0.8 are by a factor of ~ 5 and 7 larger than those 
for q 0.929 and q = I with m 2.7 Mq. Essentially 
the same results are found for the models with ALF2 



and MSI with q < 1 (see Table [IV|). By contrast, for 
the models with APR4, the total rest mass and kinetic 
energy depend weakly on the mass ratio for m = 2.7 Mq^ 
and they are always larger than those with ALF2, H4 
and MSI for m = 2.6 - 2.SMq. These facts indicate 
that (i) for relatively stiff EOS such as ALF2, H4 and 
MSI, the asymmetry of binary neutron stars enhances 
the efficiency of the angular momentum transport via the 
tidal torque and increases the total amount of the ejected 
material, and (ii) for a relatively soft EOS, APR4, which 
yields a small-radius neutron star, the total amount of 
the ejected material is always large irrespective of the 
mass ratio for the canonical total mass ~ 2.6 - 2.SMq. 
This is probably because for APR4, the shock heating in 
the early evolution stage of the HMNSs, in which they 
quasiradially oscillate with significant amplitude, plays a 
primary role in the mass ejection irrespective of the mass 
ratio; indeed, a large mass ejection is observed in the first 
1-2 ms after the onset of the merger. 

The total rest mass and kinetic energy of the ejected 



material depend also on the total mass of the system, and 
the degree of the dependence depends also on the EOS: 
For APR4, these quantities are larger for mor e m assive 
system irrespective of the mass ratio (see Table IV). This 



property is consistent with the fact that these quantities 
are larger for an EOS that yields compact neutron stars. 
Namely, for the larger mass, the system can be in general 
more compact for the binary neutron stars, and also the 
formed HMNS can be more compact, more rapidly ro- 
tate, and quasiradially oscillate with a larger amplitude. 
Thus, the mass ejection is enhanced through the angular 
momentum transport via the tidal torque and the shock 
heating. 

For H4, the similar results are obtained except for 
model H4-140140. for which the rest mass and kinetic 
energy of the ejected material is smaller than those for 
models H4-135135. The possible reason is that for H4- 
140140 (for which a black hole is formed ~ 10 ms after the 
onset of the merger), the HMNS formed is compact (i.e., 
it can trap the material in its vicinity) and moreover, its 
shape (it is not a sharp ellipsoid) could be unsuitable for 
efficiently exerting the torque to the surrounding mate- 
rial. 

For MSI for which neutron stars and HMNSs are not 
very compact, the quantities of the ejected material do 
not change very systematically. As mentioned above, for 
g < 1, the rest mass and kinetic energy increase with 
the decrease of q in the same manner as that for other 
EOSs. However, for the equal-mass case {q = 1) with 
m = 2.6 and 2.7 Mq^ the ejected rest mass and kinetic 
energy are quite large by contrast to that for q = 0.929 
and m = 2.7 Mq. This indicates that for this system, 
the HMNS is not very compact and does not trap the 
material strongly, and hence, angular momentum trans- 
port due to the torque exerted by the HMNS and shock 
heating, which are not as efficient as those in softer EOSs, 
are still large enough to overcome the trapping effect. For 
m = 2.7M0, in particular, the difference in the results of 
q = 1 and q = 0.929 is quite large. The possible reason 
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is that (i) for the equal-mass case, the ampUtude of the 
quasi-radial oscillation is by a factor of ~ 2 larger than 
that for the unequal-mass case, and thus, a larger amount 
of the materials are likely to gain the escape velocity; (ii) 
the shape of the HMNS is quite different between two 
models: For g = 1, a sharp ellipsoid is formed, and it 
appears to play a substantial role for a coherent angular 
momentum transport from the HMNS to the material 
surrounding it. By contrast, for q = 0.929, the shape is 
not a clear ellipsoid but a pear-shaped asymmetric ob- 
ject, and hence, the transport process does not appear to 
proceed efficiently. 

For MS 1-140 140, the quantities of the ejected material 
is much smaller than those for MS 1-1 35 135 and as small 
as those for MSl-130140. The possible reason is that the 
HMNS formed for MSl-140140 is slightly more compact 
than that for MSl-135135, and the mass ejection is sup- 
pressed by the trapping effect. This suggests that for 
these mass ranges, a slight change in the compactness 
significantly affects the efficiency of the mass ejection. 



B. Properties of the merger and mass ejection: 
black hole formation case 



We briefly summarize the properties of the ejected ma- 
terial for the case that a black hole is promptly formed 
after the onset of the merger. In this study, the prompt 
formation of a black hole occurs only for APR4 with the 
total mass 2.9M0. 

For these models, the mass ejection primarily proceeds 
at the instance of the merger, i.e., during a short duration 
before the formation of a black hole. Because a black hole 
is promptly formed, a region shock-heated at the collision 
of two neutron stars is soon swallowed by the black hole, 
and thus, the shock heating does not play a primary role 
in the mass ejection. A significant mass ejection occurs 
for the case that the mass asymmetry is present, and the 
mass ejection is induced primarily by a tidal torque. In 
the presence of mass asymmetry, the less-massive neutron 
star is tidally elongated during the merger, and a frac- 
tion of the tidally elongated neutron-star material gains a 
sufficient torque from the merged object just before the 
formation of a black hole and gets the escape velocity. 
For models APR4-140150 and APR4-130160, this gain of 
the angular momentum is large enough to eject materials 
of rest mass ~ 6 x lO^^M© and 2 x IO^^Mq, respectively 



(see Fig. 18). In these cases, disks are also formed, and 



their rest mass (for material bounded by the black hole) 
is O.O3M0 and O.OO2M0, respectively. The values for the 
mass ejection depend only very weakly on the grid res- 
olution with the fluctuation within 10 - 20% level (see 
Appendix A). The reason is that strong shocks do not 
play an important role in the mass ejection mechanism. 

The average velocity of the ejected material for these 
cases is ^ 0.3c and larger than that in the case of the 
HMNS formation. The reason is that the mass ejection 
is caused primarily by the tidal interaction at the on- 



set of the merger, and for this case, the induced velocity 
is larger than that by subsequent longterm shock heat- 
ing. Because the tidal interaction plays a primary role, 
the material is ejected primarily in the direction of the 
equatorial plane. The motion to the z direction is also 
induced by shock heating that occurs when spiral arms 
surrounding the black hole collide each other. However, 
this is a secondary effect. Hence, for the case that a black 
hole is promptly formed from an asymmetric binary, the 
value of ^0 is 30 - 35° which is smaller than those for the 
case of the HMNS formation for which 6>o = 40 - 50°. 
For the equal-mass binary, the total rest mass of the 



ejected material is quite small ^ 10 ^M© (see Fig. |18[), 



because of the absence of the asymmetry and of the lack 
of the time during which the material located in the outer 
region gains the torque from the merged object (note that 
most of the fluid elements of binary neutron stars just 
before the onset of the merger do not have the specific 
angular momentum large enough to escape from the black 
hole formed [21]). In this case, the disk mass surrounding 
the black hole is also quite small, ~ 10~'^Mq. This is 
consistent with our previous finding [25 . 

Figure 19 plots the gravitational waveforms for APR4- 
130160 and APR4-140150. For these models (also for 
APR4-145145), the gravitational waveform is character- 
ized by the inspiral waveform and subsequent ringdown 
waveform. The frequency of gravitational waves mono- 
tonically increases and eventually reaches the value of the 
fundamental quasinormal mode of the formed black hole. 
For all three cases, the frequency of gravitational waves 
associated with the quasinormal mode is 6.55 ±0.05 kHz 
(the mass and spin of the formed black holes are 2.8 - 
2.83M0 and 0.77 - 0.78, respectively), which agrees with 
the frequency of the quasinormal mode analytically de- 
rived [52 . Because HMNS is not formed, no feature for 
the quasiperiodic oscillation associated with the HMNS 
formation is found. 



V. SUMMARY AND DISCUSSION 

A. Summary 

We reported our latest numerical-relativity studies for 
the material ejected in the merger of binary neutron 
stars. We explored the properties of the ejected material 
for a variety of EOSs, total masses, and mass ratios of bi- 
nary neutron stars, and found the following facts. First, 
we summarize the results for the case that a HMNS is 
formed: 

• For the canonical total mass of the binary neutron 
stars 2.6 - 2.8M0, the total rest mass and kinetic 
energy of the ejected material are approximately in 
the range 10""^ - lO-^M© and 10"^^ - 10^^ ergs, re- 
spectively, depending on the EOS, total mass, and 
mass ratio of binary neutron stars. 

• The total rest mass and kinetic energy of the 
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FIG. 17: The same as Fig. [3] but for high-mass and unequal-mass model APR4-130160. The filled black circles in the middle 
and right panels of the top row denote black holes. 
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FIG. 18: M*esc and T*esc as functions of t — tAH for mod- 
els APR4-130160, APR4-140150, APR4-145145. Here, Uh 
denotes the time at which an apparent horizon is formed. 



ejected material depend strongly on the EOS. They 
are in general larger for binaries composed of an 
EOS that yields compact (small-radius) neutron 
stars (e.g., APR4). They also depend on the ef- 
ficiency of shock heating (i.e., depend on Fth)^ For 
many cases, smaller values of Fth results in a larger 
ejected mass and kinetic energy. 

• The total rest mass and kinetic energy could de- 
pend also strongly on the mass ratio of binary neu- 
tron stars. The dependence on the mass ratio is in 
particular strong when the radius of neutron stars 
is relatively large (i.e., for stiff EOSs such as H4 
and MSI). For many models in such EOSs, the 
ejected rest mass and kinetic energy are larger for 
smaller mass ratios (for more asymmetric binary 
neutron stars). By contrast, when the EOS is rela- 
tively soft (i.e., for APR4), the dependence of these 
quantities on the mass ratio is weak. 

• The total rest mass and kinetic energy depend also 
on the total mass of binary neutron stars. For many 
models in the present study, these quantities are 
larger for the larger total mass irrespective of the 
EOS. However, there are some exceptions for a class 
of EOS and for a massive system. 

• The total rest mass ejected is in the wide range 
between - 5 x IQ-"^ (H4) and - 7 x IO'^Mq 
(APR4) for equal-mass binaries with the total mass 
m = 2.7 Mq. For the unequal-mass case with 
q ~ 0.8, it is in a rather narrow range between 
- 3 X 10-3 (MSI and H4) and - 8 x IO'^Mq 
(APR4) for m = 2.7 Mq. This indicates that the 
ejection of the material is induced by two different 
processes, the torque exerted the HMNS and the 
shock heating. 

• The total kinetic energy is also in the wide range 
between - 10"^^ ergs (MSI and H4) and - 10^^ ergs 



(APR4). The typical maximum velocity of the 
ejected material is 0.5 - 0.8c, and the average veloc- 
ity is 0.15 - 0.25c. For EOSs that yield a compact 
neutron star, the velocity of escaping material is 
larger. 

In our present study, a black hole is promptly formed 
for three models; APR4-145145, APR4-140150, and 
APR4-130160. For the case of the prompt black-hole 
formation, a significant mass ejection occurs only for the 
asymmetric binaries. For APR4-145145, the rest mass 
of the ejected material is ~ 10~^Mq. By contrast, it 
is rather large for APR4-140150 and APR4-130160 as 
6 X 10""^ M0 and 2 x 10'^ Mq. Thus, a massive and 
asymmetric binary can eject a large amount of the mate- 
rial even if a black hole is promptly formed (even in the 
absence of a HMNS). For this case, the average velocity of 
the ejected material is larger than that for the case of the 
HMNS formation, and thus, a unique feature may be seen 
in the observed electromagnetic signal (see next section) . 
However, we should note that the prompt black-hole for- 
mation will occur only for a large total mass, because the 
latest discovery of a high-mass neutron star PSR J 1614- 
2230 with mass 1.97 ± O.O4M0 [29] indicates that the 
EOS should be rather stiff and with such stiff EOSs, a 
HMNS is the canonical outcome for the canonical-mass 
binary with m = 2.6 - 2.8M0. 

The distribution of the matter around the remnant ob- 
ject depends strongly on the merger process. For the case 
that a HMNS is formed, a dense (physical) atmosphere 
(including the ejected material) is formed around it (see 
Figs. |3]-[5|. The atmosphere is distributed in a weakly 
anisotropic manner, and thus, the matter is present even 
in the spin axis of the HMNS. Such dense atmosphere 
will be present even after the HMNS collapses to a black 
hole. HMNS or black hole subsequently formed will emit 
a huge amount of neutrinos (e.g., [27 ), and may drive a 
jet from the central region through the fireball produc- 
tion via neutrino-antineutrino pair annihilation. To drive 
a SGRB for this case, however, the jet has to penetrate 
the dense atmosphere and ejected material. Whether it 
is possible or not is a question to be answered by the 
future research. By contrast, for the case of the black 
hole formation, the atmosphere is not very dense around 



the spin axis of the black hole (see Fig. 17). Thus, for 
this case, a SGRB would be driven, if an energetic jet is 
launched as studied, e.g., in [53 . 



B. Discussion 

We here briefly argue possible electromagnetic signals 
emitted by the material ejected from the merger of bi- 
nary neutron stars, referring to the numerical results in 
the present work. As already mentioned, a recent discov- 
ery of a high-mass neutron star PSR J1614-2230 suggests 
that the maximum mass of spherical neutron stars should 
be larger than 1.97 ± O.O4M0 [29]. This indicates that 
a long-lived HMNS would be a canonical outcome of the 
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FIG. 19: The same as Fig.[9]but for models APR4-130160 (left) and APR4-140150 (right). 



merger of binary neutron stars, if the binaries were com- 
posed of neutron stars of canonical mass of 1.3 - IAMq 
with the total mass ~ 2.6 - 2.SMq [26^. The present nu- 
merical results indicate that from the long-lived HMNS, 
a fraction of the material could be ejected with large ki- 
netic energy. 

References [T]-[3l [18] discuss the signals by the radioac- 
tive decay of r-process nuclei, which would be produced 
from the neutron-rich material in the outflow [T2]-[17], 
and subsequently decay and emit a signal that may be ob- 
servable by current and future-planned optical telescopes 
such as PTE [54 , Pan-STARRs [55], and LSST [56]. In 
this scenario, the typical duration of the peak luminosity 
is of order a day or less as [1 



t 



peak 



: O.ld 



0.2 



-1/2 



V 10-3 

and the associated peak luminosity is 

/eff 



1/2 



(29) 



i^peak ~ 7 X 10"^^ ergs/s 



0.2 



3 X 10-6 

1/2 



10-3Mc 







1/2 



(30) 



where /eff denotes the conversion rate of the energy per 
rest-mass energy in the ejected material through the ra- 
dioactive decay process, which is ~ 3 x 10"^ according 
to the results of [3]. Pqc is the typical velocity of the 
ejected material. The result of [3 suggests that if the 
total ejected mass is ^ lO^^M©, the signal will be de- 
tected by large optical surveys such as LSST for a typical 
distance to sources ~ 100 Mpc. Our numerical results in- 
dicate that /3o = 0.15 - 0.25, and the total ejected mass is 
~ 10~3 - 1O~^M0 for binaries composed of neutron stars 
with a small radius ^ 11 - 12.5 km, and ~ 0.3 x 10"^ - 
5 X 10~^Mq for binaries composed of neutron stars with 
a larger radius ~ 13.5 - 14.5 km for plausible values of 
Tth = 1.6 - 2.0. Thus, if the EOS is a rather "soft" one 
that yields a small-radius neutron star, an observable op- 
tical signal due to the radioactive decay can be expected 



with a duration of several hours. Taking into account the 
short duration of the signal, rapid follow-up searches and 
an efficient coverage for the error circle of the direction 
of the gravitational- wave events are required. If the EOS 
is a "stiff" one that yields a large-radius neutron star, 
the strength of the signal will be weaker and the dura- 
tion shorter, although it would be still possible to detect 
the signal in particular for the merger of unequal-mass 
(sufficiently asymmetric) neutron stars. 

There is also another possible channel for the electro- 
magnetic emission. According to recent studies [H [181 
[20] . the ejected material, which is in the free expan- 
sion, will sweep up the interstellar matter and form blast 
waves. During this process turning on, the shocked ma- 
terial could generate magnetic fields and accelerate par- 
ticles that emit synchrotron radiation, for a hypothetical 
amplification of the electromagnetic field and a hypo- 
thetical electron injection. The emission will peak when 
the total swept-up mass approaches the ejected mass, be- 
cause the blast waves are decelerated and transit to the 
phase in which the motion of the material is described 
by the (non-relativistic) Sedov- Taylor's self-similar solu- 
tion. The predicted duration for the synchrotron radi- 
ation depends on the total energy Eq and speed of the 
ejected material Pqc as well as the number density of the 
interstellar matter no- The duration to reach the peak 
luminosity is estimated in [4] as 



Tradio ^ 4 yrs 



En 



IQ50 gj.gg 



1/3 / X -1/3 / o \ -5/3 



no 



1 cm 



V 0.2 



(31) 

By the synchrotron radiation, a radio signal could be 
emitted as in the late phase of supernovae and the after- 
glow of gamma-ray bursts [4 . Our numerical results in- 
dicate that the typical velocity of the ejecta is Po = 0.15 
- 0.25 irrespective of the EOS and masses of neutron 
stars in binaries. However, Eq is in a wide range be- 
tween ^ 10^^ ergs and 10^^ ergs, depending strongly on 
the EOS, mass ratio, and total mass of the binaries, and 
its value is highly uncertain. Thus the predicted value 
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of '^radio IS in a wide range ~ 1 - 10 yrs, even for an 
optimistic value of no = 1 cm~^. For smaller values of 
no which is likely when the merger occurs outside the 
galactic plane, the value of Tradio is much longer. 

For the typical value of the ejecta velocity f3o ~ 0.2, 
the peak flux for the observed frequency is obtained at 
the deceleration time described in Eq. (31). Specifically, 
the peak flux may be obtained at the self-absorption 
frequency, ^ 1-2 hundreds MHz, and the typical syn- 
chrotron frequency is sub-MHz. The peak flux for a given 
observed radio-band frequency i^obs is 



^ 90 /iJy 



^0 



(r 



no 



2Q50 g^gg y v 1 cm 



D 



0.9 



0.2 



200 Mpc 



Vl.4 GRz) 



-0.75 



(32) 



where we assumed the power-law distribution of the elec- 
tron's Lorentz factor with the power 2.5. Equation (32) 



is applicable as long as the observed frequency is higher 
than the typical synchrotron and self-absorption fre- 
quency at the deceleration time, Tradio- Equation (32) 



indicates that for a hypothetical event at a distance of 
200 Mpc, Eq ^ 10^^ ergs with no = 1 cm"^ is strong 
enough to be observed by future-planned radio instru- 
ments (such as EVLA [57], ASKAP [58], MeerKAT [59], 
and Apertif for which the root-mean square value of the 
background noise for one hour observation is smaller than 
50 /iJy as shown in [4 ). Therefore, the mass-ejection 
mechanism could supply a large amount of the kinetic 
energy which generates an observable strong radio signal, 
if the EOS is rather soft (i.e., the neutron-star radius is 
fairly small) or the binary is significantly asymmetric. 

In this scenario, the duration to reach the peak lu- 
minosity and the strength of the radio signal depend 
strongly on the value of no- In nature, the value of no will 
depend strongly on the site where the merger of binary 
neutron stars happens. If it is in a galactic disk, no would 
be typically ~ 1 cm~^, whereas if it is outside a galaxy, 
the value is much smaller as ~ 10~^ cm~^. Equation 
( [32I shows that for a smaller value of no <C 1 cm~^, 
Fjy ^ 1/iJy even for Eq = 10^^ ergs. Our numer- 
ical simulation shows that the maximum value of £^0 
is at most 10^^ ergs. Therefore, for the low value of 
no ~ 10~^ cm~^, this type of electromagnetic signals may 
not be observable as a counterpart of the gravitational- 
wave signal [18 . 

We here note the following point. We used the total 
kinetic energy and average velocity as £^0 and /3o, when 
estimating the radio flux estimated in Eq. (32). How- 



ever, as we found in this paper, the ejected material has 
a wide range of the velocity, and the amount of the ki- 
netic energy which the material of a given value of has 
depends on the value of Pq. Thus, each material has a 
different deceleration time and flux F^, and therefore, the 
light curve will have a complicated structure depending 
on the distribution of the differential mass as a function 
of the ejecta velocity (see also [19 ). In a subsequent pa- 



per, we plan to study the luminosity curve in more detail 
following [20 . 

V A[ the properties of the 



As summarized in Sec. 



ejected material depend strongly on the EOS, mass ra- 
tio, and total mass of the binary. This suggests that 
the observed electromagnetic signal depends on them. 
The observation of gravitational waves in the inspiral 
phase of binaries, which will be observed in the advanced 
gravitational- wave detectors [5 , will carry the informa- 
tion of the mass ratio and total mass. The observation 
of gravitational waves from the final inspiral phase and 
HMNS could constrain the EOS of neutron stars. To- 
gether with these information by the gravitational-wave 
observation, the observation of the electromagnetic sig- 
nals will be used for clarifying the dynamics of the bi- 
nary merger and ejected material. In addition, clarifying 
the spectroscopic properties of electromagnetic emission 
associated with the decay of r-process nuclei may be in- 
teresting. If there is a characteristic emission/absorption 
feature in this emission, it will be helpful to determine 
the cosmological redshift of the source event through a 
spectroscopic observation. The redshift determined will 
be subsequently used for determining the distance to the 
source (and thus, Hubble constant) and the physical mass 
of binary systems through the analysis of observed inspi- 
ral gravitational waves [60 . However, for the clarifica- 
tion, detailed theoretical studies for the ejected material 
and electromagnetic radiation are necessary for a vari- 
ety of the EOS, mass ratio, and total mass. These are 
new and interesting tasks for the community of numerical 
relativity. 
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Appendix A: Convergence 

The rest mass and kinetic energy of the ejected ma- 
terial with different grid resolutions for selected models 
are listed in Table [V| This shows that for the unequal- 
mass models, a convergence is well achieved (even for 
mi = I.3M0 and m2 = I.4M0), whereas for the equal- 
mass models, the convergence is poor. In this case, the 
results with > 60 could be modified by a factor of ~ 2 
from the results with N = 60. As mentioned in Sec. IV A[ 
the convergence is poor for the case that a strong shock 
is formed at the merger and it plays a primary role in the 
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mass ejection. The possible reason for this is that (i) the 
shocks are always computed by the first-order accuracy, 
and hence, the accuracy is low and (ii) the ejected mass 
is a tiny part of the entire system, and hence, a random 
error for the entire system computed with a low accuracy 
significantly (and randomly) affects a tiny amount of the 
ejected material. By contrast, when the tidal torque and 
hydrodynamical torque exerted by the HMNS play an 
important role in the mass ejection, the convergence is 



good. 

The averaged frequency of gravitational waves emitted 
by HMNSs is also listed for three grid resolutions. This 
shows that the frequency is obtained within the error of ^ 
0.1 kHz, which is smaller than the physical dispersion of 
the frequency associated with the quasiradial oscillation 
of the HMNSs and the secular change of the density and 
velocity profiles of the HMNSs caused by the angular 
momentum transport. 
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